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ABSTRACT 

We report radio polarization observations of G319.9-0.7 (MSC 319.9-0.7) at 3 and 6 cm obtained 
with the Australia Telescope Compact Array. The source shows a highly elongated morphology 
with the energetic pulsar J1509-5850 located at the tip. We found a flat radio spectrum of index 
a = —0.26 ± 0.04 and a high degree of linear polarization. These results confirm G319.9-0.7 as a 
bow-shock pulsar wind nebula. The polarization maps suggest a helical magnetic field trailing the 
pulsar, with the symmetry axis parallel to the system's inferred direction of motion. This is the 
first time such a field geometry has been seen in a bow-shock nebula, and it may be the result of an 
alignment between the pulsar spin axis and its space velocity. Compared to other bow-shock examples, 
G319.9-0.7 exhibits very different properties in the field structure and surface brightness distribution, 
illustrating the large diversity of the population. 

Subject headings: pulsars: individual (PSR J1509-5850) — stars: neutron — stars: winds, outflows 
- radio continuum: ISM — ISM: individual objects (MSC 319.9-0.7) 



1. INTRODUCTION 

Pulsars lose a significant fraction of their rotational 
energy through their rclativistic winds. The consequent 
interactions with the ambient medium result in broad- 
band synchrotron emission. These structures are col- 
lectively referred to as pulsar wind nebulae (PWNe). 
The properties of a PWN depend strongly on its evo- 
lutionary state and environment. Since a pulsar typi- 
cally has a space velocity of a few hundred kilometers 
per second, it will eventually escape its natal supernova 
remnant (SNR) and travel supersonically in the inter- 
stellar medium (ISM). In such cases, the pulsar outflow 
can be confined by ram pressure, resulting in a bow- 
shock nebula. Multiwavelength observations have iden- 
tified several bow-shoc k PWN systems (see review by 
IGaensler fc Sland 12006). The best-studied example is 
'the Mouse' (G359.23-0.82), in which the X-ray and radio 
emissions can be well-modeled by a bright head coinci- 
dent with the pulsar, a 'tongue' region corresponding to 
the wind termination shock, and an elon gated tail asso- 
ciate d with the post shock flow material (Gaenslcr et al. 
[200l . 

Studies of the Mouse and other bow-shock PWNe 
have set the stage for su bsequent theoreti c al an d nu- 
merical modeling efforts. iRomanova et al.l (|2005[ ) con- 
sidered an axisymmetric case of a pulsar traveling along 
its spin axis direction. They proposed an azimuthal B- 
field geometry for bow-shocks, and derived the analyti- 
cal shape of a pulsar tail. With a similar aligned con- 
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figuration, iBucciantini et ahl ([20051 ) carried out the first 
relativistic magnetohydrydynamic (MHD) simulations of 
bow-shocks to study the e f fects o f different wind mag- 
netization. I Vigelius et al.1 (|2007t ) relaxed the assump- 
tion of alignment and presented three-dimensional non- 
relativistic hydrodynamic simulations to illustrate the 
dependence of bow-shock morphology on pulsar orien- 
tation and ISM gradient. 

Observationally, PWNe at radio frequencies are char- 
acterized by flat spectra with spectral indes0 a ~ —0.3 
to 0. and high (> 1 0%) degrees of linear polariza- 
tion ([Kaspi et al.1 12006) . Due to the long synchrotron 
cooling time of the radio-emitting particles, a radio 
PWN can act as a direct calorimeter to reflect the sys- 
tem's integrated history over a long period of time. 
More importantly, radio polarimetry provides a powerful 
probe of the magnetic field geometry of a PWN, which 
is relatively less studied. High resolution polarization 
measurements have only been carried out on a hand- 
ful of bow-shock s y stems, e.g. IC 443 and the M ouse 
(|Qlbert et alJl200lt lYusef-Zadeh fc Gaenslerl 120051) . In 
particular, VLA observations of the Mouse revealed a 
highly ordered field wrapping around the bow shock at 
the apex and extending parallel to the ne bular axis down- 
stream (jYusef-Zadeh fc Gaensler] [2005h. in contrast to 
the theoretical prediction by IRomanova et al.l (|2005| ). 

To generalize our understanding of PWNe, we need 
to study more examples. In this paper, we report radio 
observations of a bow-shock nebula powered by the en- 
ergetic pulsar J1509-5850. This pulsar was discovered 
in th e Parkes Multibeam Pulsar Survey (jKramer et al.1 
l2003h . It has a spin period P = 88.9 ms, a high spin- 
down luminosity E = 5.1 x 10 35 ergss -1 and a relatively 

4 The spectral index a is related to the flux density Sv and 
frequency v by S v oc v a . 
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young characteristic age t c = P/(2P) = 1.5 x 10 5 yr. 
As with other energetic pulsars, it has recently been de- 
tected in 7-rays by the Fermi Gamma-ray Space Tele- 
scope (|Weltevreddl20Tl . The pulsar dispersion measure 
(DM) of 137.7pccm -3 suggests a distance ranging from 
2.6 to 3. 8 kpc according to different Galactic free electro n 
models (|Tavlor fc Cordesl [Tool ICordes fc Laziol 12^ . 
Throughout this work, we will adopt a pulsar distance 
d = 3 c?3 kpc. Some ~ 4' southwest of the pulsar, there 
is an elongated radio source MSC 319.9-0.7 (hereafter 
G319. 9-0.7) which was identified by iWhiteoak fc Greenl 
(1996) as a candidate SNR based on observations with 
the Molonglo Observatory Synthesis Telescope (MOST). 
The radio source has a size of 7' x 2', and consists of a 
brig ht central bulge an d a clump in the south. 

iHui fc Beckerl (|2007l ) reported an observation of the 
field of PSR J1509-5850 with the Chandra X-ray Ob- 
servatory. The X-ray image revealed a tail-like PWN 
extending from the pulsar in the same orientation as 
G319.9-0.7, suggesting that it is a bow-shock nebula. 
At 3 kpc, the tail has a physical scale ~ 5pc, one of the 
longest X-ray tails ever observed. The PWN has an X- 
ray luminosity Lx ~ 10 33 ergs s -1 between 0.5 and 8 keV, 
corresponding to a high efficiency -q = Lx / E w 2 x 10 -3 . 
IHui fe Beckerl (|2007f l argued that G319.9-0.7 is too small 
to be a SNR associated with the pulsar, and concluded 
that i t is most likely a backgr ound obje ct. O n the other 
hand, iKargaltsev et al.l (|2008l . hereafter IK08D presented 
a detailed study using the same Chandra data set, and 
suggested that G319.9-0.7 could be the radio counterpart 
of the X-ray PWN. They found an average flow speed 
> 5000 km s -1 and an equipartition field of a few times 
10 -5 G in the X-ray-emitting region. 

Motivated by these previous studies, we have carried 
out high resolution radio imaging using the Australia 
Telescope Compact Array (ATCA) to investigate the na- 
ture of G319.9-0.7. Observations and data reduction are 
described in We present the results in fJ3]and discuss 
the physical implications in £|4] Finally, we summarize 
our findings in fj5] 

2. OBSERVATIONS AND DATA REDUCTION 

ATCA imaging of G319.9-0.7 was carried out in two 
12 hr observations with the EW 367 and 750C array con- 
figurations on 2007 Nov 16 and Dec 1, respectively. To 
obtain a uniform u-v sampling and to boost the signal- 
to noise ratio, the longest baselines were excluded in our 
analysis, giving array spacings from 46 m to 750 m. The 
observations were made at 3 and 6 cm (center frequen- 
cies at 8.64 and 4.8 GHz, respectively) simultaneously 
with all Stokes parameters recorded. At each wave- 
length, the total bandwidth of 128 MHz was divided into 
32 x 4 MHz overlapping spectral channels. After dis- 
carding the edge channels and channels affected by self- 
interference, and merging the overlapping adjacent chan- 
nels, a usable bandwidth of 104 MHz is left, which is split 
into 13 x 8 MHz channels. The entire field was covered 
by a 9-pointing mosaic patten with pointings spaced by 
3', smaller than the FWHM of the primary beam even 
at 3 cm. 

The flux density scale is set by observations of the pri- 
mary calibrator, PKS B1934-638. In addition, a sec- 
ondary calibrator, PKS B1511-55, was observed every 
30 minutes to determine the antenna gains. All our data 



reduction was performe d with standard te chniques using 
the MIRIAD package (jSault et all 119951 ). We first ex- 
amined the data carefully to reject outlier visibility data 
points and scans during poor atmospheric phase stabil- 
ity. We then determined gain, bandpass, flux and po- 
larization calibrations for the observations, and applied 
these to the source visibilities. Mosaiced images in each 
polarization were formed using the mult ifrequency syn- 
thesis technique (jSault fc Wieringal [l~994') to improve the 
u-v coverage and with uniform weighting to minimize 
the sidelobes. Images in Stokes I, Q and U were decon- 
volved simultane ously using a ma ximum entropy algo- 
rithm (PM0SMEM; iSault et all 119991 ). and restored with a 
Gaussian beam of FWHM 17"xl5", corresponding to the 
diffraction limit at 6 cm with the 750 m baseline. Finally, 
we convolved the 3 cm map to identical spatial resolution 
as the 6 cm one. Our final maps in Stokes 7, Q and U 
at 3 cm have rms noise of 0.2 mJy beam -1 , very close to 
the theoretical level of 0.18 mJy beam -1 . At 6 cm, the 
noise level is 0.15 mJy beam -1 in the Stokes I map, and 
0.1 mJy beam -1 in the Stokes Q and U images, again 
consistent with the theoretical value of 0.1 mJy beam -1 . 

Using the Stokes Q and U images, we generated maps 
of the polarized intensity and position angle (PA) with 
the t ask IMPOL, which includes corrections for Ricean 
bias (jWardlc fc Kronberelll974D . To obtain the intrin- 
sic PAs of the polarization vectors, we reprocessed the 
20 cm ATCA data from the Souther n Galactic Plane 
burvey (SGPS; lHaverkorn eFal1l2"00"l on a channel-by- 
channel basis, applied the rotation measure (RM) synthe- 
sis technique (|B rent i ens fc de Bruvnl I2005T ) to map the 
foreground Faraday rotation, and then derotated the 3 
and 6 cm polarization vectors accordingly. 

In addition to the radio data, we have also processed 
the 40 ks Chandra X-ray observ ation (Ob sID 3513) re- 
ported by IHui fc Beckerl (|2007f ) and IK08I for a multi- 
wavelength comparison. After the standard pipeline pro- 
cessing, we generated an exposure-corrected image in the 
0.5-7 keV energy range and applied a smoothing to 6" 
resolution to maximize the signal-to-noise ratio. 

3. RESULTS 

3.1. Nebular Morphology 

Radio continuum images of G319.9-0.7 at 3 and 6 cm 
are presented in Figure [T] The source exhibits a very sim- 
ilar morphology at both wavelengths; it is highly elon- 
gated and extends over 11' with a clear symmetry axis 
along the PA 200° (measured from north through east). 
We note that since the shortest array spacing in the data 
is only 46 m, the sensitivity of the 3 cm intensity map is 
expected to drop for angular scale beyond ~ 3', resulting 
in obvious sidelobes in FigureQJb). PSR J1509-5850 was 
not detected in our observations, but its position, which 
is marked by the cross in the figure, coincides with the tip 
of the nebula. At the pulsar location, the radio emission 
is very faint. It gradually brightens towards the south 
and shows a cone-like morphology with an opening angle 
~ 35°. The emission peaks at a central bulge 4f5 from 
the pulsar and attains a maximum width of 3'. Beyond 
that, the nebula fades and narrows. At 8' southwest of 
the pulsar, there is a bright, unresolve d radi o source with 
a faint X-ray counterpart reported by IK081 which could 
possibly be a background radio galaxy. Further south, 
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the nebula changes slightly in orientation (~ 10°) and 
narrows to a collimated faint tail. Finally, beyond ~ 11' 
from the pulsar, the nebular surface brightness drops to 
near the noise level of the radio maps, although there is 
some hint of extension for another 1-2 arcminutes. 

Figure [2] shows the exposure-corrected Chandra X-ray 
image, overlaid with radio contours from the 6 cm ra- 
dio intensity map in Figure QJa). The X-ray and ra- 
dio emissions near the pulsar are well-aligned and are 
bounded by a common envelope. We followed IK08I by 
modeling the envelope with a parabola and obtained 

1 /2 

x = 3" [{z + 0. 5)/l J , where x and z are the angular- 
distance (in arcseconds) from the pulsar, measured per- 
pendicular and parallel to the nebular symmetry axis, re- 
spectively. As a compariso n, ou r result is slightly wider 
than the one reported by IK081 which may due to the 
lower resolution of the radio map than the X-ray im- 
age. Figure [2] also indicates an anti-correlation between 
the radio and X-ray surface brightness. While the for- 
mer increases with distance from the pulsar, the latter 
peaks at the pulsar position and fades to the south. At 
1.'5 downstream from the pulsar, the X-ray tail suddenly 
narrows by a factor of 2. We found some hint of a radio 
ridge, as indicated by the middle contour in Figure [21 
which may correspond to the narrowed X-ray tail. Fur- 
ther south, the X-ray emission becomes more diffuse and 
seems to extend to the radio peak. Deeper observations 
are needed to confirm this. Figure [3] shows the surface 
brightness profiles of the radio and X-ray emissions, il- 
lustrating the features described above. 

3.2. Polarization Properties 

After correction for the foreground Faraday rotation 
(see below), the strength and intrinsic projected orien- 
tation of the polarization B-vectors of G319.9-0.7 are 
plotted in Figure 21 revealing a highly ordered magnetic 
field structure. The B- vectors in the north runs along the 
northwest-southeast direction with a PA of 120°, nearly 
perpendicular to the symmetry axis of the nebula. South 
of the bulge, the vectors switch direction abruptly, and 
show a good alignment with the nebular axis. 

The fractional linear polarization map of G319.9-0.7 
at 6 cm is shown in Figure [SJ The source is highly lin- 
early polarized and the degree of polarization appears to 
be systematically higher around the edges than in the 
interior. For instance, it is over 40% polarized along the 
southern edge and > 30% in the north, but only < 25% 
at the central bulge. The map at 3 cm, which is not 
shown here due to low signal-to-noise ratio, suggests a 
very similar pattern. 

Figure [H] presents the results from RM synthesis of the 
20 cm data, indicating a small magnitude of RM for most 
regions. The central bulge has |RM| < 20radm~ 2 , and 
the RM slightly increases to +50radm -2 towards the tip 
in the north. The measurement uncertainties in the map 
are of a similar order (< 50radm~ 2 ). Since the compact 
source in the south is very weakly polarized, we cannot 
determine if it has a different RM than the rest of the 
nebula. Using this RM map, we derotated the 3 and 6 cm 
polarization vectors to their intrinsic orientation, which 
are plotted in Figure 2J Due to the small RM values, the 
corrections at both wavelengths are generally small, e.g. 
a RM of 50radm~ 2 would correspond to a derotation of 



3° at 3 cm and 11° at 6 cm. As a note, the measurement 
errors in the PA of the 3 and 6 cm vectors are about 8° 
and 4°, respectively. 

3.3. Radio Spectrum 

After subtracting the unresolved source in the south, 
G319.9-0.7 has flux densities of 0.29 ± 0.02, 0.43 ± 0.02, 
0.60 ± 0.05 and 0.69 ± 0.05 Jy at 3, 6, 20 and 36 cm, 
respectively, as obtained from our data (3 and 6 cm), 
SGPS (20 cm: lHaverkorn et al.l l200ll and MOST (36 cm; 
iWhiteoak fc Greenlll996f ). The resulting radio spectrum 
is plotted in Figure[7{a). With the shortest array spacing 
of 46 m, the 3 cm intensity map is only sensitive to angu- 
lar scales smaller than 2. '6. Therefore, the flux measure- 
ment at 3 cm is likely underestimated and excluded this 
data point to deduce a spectral index a = —0.26 ± 0.04 
from 3 other bands. The compact source has a substan- 
tially steeper radio spectrum of a = —0.96 ± 0.02, sug- 
gesting th at it could be an unr e lated background source. 
Following iGaensler fc Wallace! (|2003i ). we resampled the 
6 cm data and co mpared to the 3 cm map thro ugh spec- 
tral tomography (jKatz-Stone fc RudnickllT997T ). The re- 
sult indicates no significant spatial variation in the spec- 
tral index along the nebula. Fi nally , a comparison to 
the X-ray spectrum reported by IK08I is plotted in Fig- 
ure[7Jb). The broadband spectrum suggests a break near 
~ 10 13 Hz with a change in spectral index Aa « 1. In- 
tegrating the radio spectrum below the break frequency 
gives a radio luminosity Lr — 8.7 x 10 33 ergss -1 , nearly 
an order higher than the X-ray luminosity. 

4. DISCUSSION 
4.1. Nature of G3 19. 9-0. 7 

Our results clearly show that G319.9-0.7 is a non- 
thermal radio source with a flat spectrum (a ~ —0.3) and 
a high degree of linear polarization (30-40%), thereby es- 
tablishing the PWN nature of the source. Moreover, the 
location of G319.9-0.7 respect to PSR J1509-5850 and 
its alignment with the X-ray tail strongly suggest that 
it is a radio counterpart of the bow-shock nebula. Pre- 
vious Chandra X-ray observation in dicat es a projected 
bow-shock standoff distance ~ 0'.'5 (K08), which is too 
small to be resolved by our radio maps. Using this valu e, 
the theoretical shape of a bow-shock (e.g. I Wilkirill 1 996H is 
about 5 times narrower than the parabola we obtained. 
This is likely the result of deceleration in the flow down- 
stream (see £ )4. 21 below). 

Due to long synchrotron cooling timescale, the radio 
emission of a PWN traces the pulsar motion over a large 
distance, reflecting the path of the system. Therefore, 
in the following discussion, we call G319.9-0.7 a radio 
'trail' rather than a 'tail', and its orientation implies a 
pulsar proper motion direction towards the northeast at 
PA ~ 20°. At the pulsar distance 3 d$ kpc, the radio trail 
has a physical exten t over 10 ^3 pc. For a space velocity 
~ 300 km s -1 (K08), if the pulsar is nearly as old as its 
spin-down age t c — 1.5 x 10 5 yr, then it would have trav- 
eled a distance 8 times longer than the radio trail. Even 
accounted for the projection, we think it is unlikely that 
the radio PWN reflects the entire history of the system, 
and the pulsar birth site could be much further south. 
At 4° southwest of the pul sar, there is an op en cluster 
Ruprecht 112 fC1453-623: lAlter et all WWm lining up 
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exactly with the nebular axis. However, this cluster is 
poorly studied and the distance is unknown. If this is 
the birth site of the pulsar, the large angular separation 
would require a high pulsar velocity of 1400 km s -1 , an 
older age, or a closer source distance. A direct measure- 
ment of the pulsar proper motion could provide more 
insights into the birth site. 

The small RM towards the source is not unexpected, 
since the line of sight to the source is near the inter- 
arm region between the Crux and Norma Galactic spiral 
arms, where the RM contributed by the ISM is small 
(|Han et al.l 120061 ). As their work suggests, our result is 
consistent with a source distance below 4kpc. The ra- 
tio between the PWN RM (< 20radm~ 2 ) a nd the pul- 
sar DM (137. 7pc cm" 3 ; iKramer et al.l I2003h indicates a 
small mean field along the line of sight 



< B» > = 1.232 — — uG < 0.2 uG 

ll DM p ~ p 



(1) 



Finally, we note that the compact sour ce in the south 
has a flux density of O.Uy at 20 cm. Hop kins et al.l 
(2003) found that radio sources with at least this flux 
level have a number density of 2.0 per square degree in 
the sky. Therefore, within the PWN area ~ 0.01 deg 2 , 
the chance probability of having such a background 
source is about 2%, not completely negligible. 

4.2. Magnetic Field Strength and Flow Structure of 
G3 19. 9-0. 7 

Figure [Tfb) shows that the X-ray synchrotron radia- 
tion of the PWN could extend down to 10 13 Hz without 
a spectral break. According to K0|| this value implies 
an equipartition held strength ~ 100 fiG near the pulsar 
and a projected flow speed ^flow.x ~4x 10 5 krns" 1 . To 
estimate the magnetic field downstream beyond the X- 
ray tail, we consider the radio emission in the southern 
half of the nebula beyond 4f5 from the pulsar. We inte- 
grate the spectrum in the range 10 7 — 10 13 Hz to obtain 
a luminosity L = 4.4 x 10 33 d§ ergss -1 (after excluding 
the background source) , with the corresponding emission 
volume^ = 2.0xl0 57 d| cm . This gives an equipartition 
field 

B cq = [Qirc 12 {l + ^L^V- 1 ]^ G w20d^(l+fc)7$-? 

(2) 

where k is the ion to electron energy ratio, < $ < 1 
is the volume filling factor, and c\2 is a constant weakly 
depends on the spectral index ()Pacholczvk| [T970)Fl As- 
suming k — and $ = 1, we obtained B cq = 20 /iG, 
a few times lower than that of X-ray-emitting region. 
We note that this result is independent of the lower fre- 
quency limit, and insensitive to the upper limit, e.g. the 
range 10 7 - 10 11 Hz yields B cq = 15 /iG. 

While the flow speed in the radio PWN cannot be de- 
termined precisely without knowing the maximum en- 
ergy of the synchrotron radiation, the absence of X-ray 
emission can still provide some rough estimates. From 
the point where the X-ray tail vanishes, the radio nebula 
extends 7' further downstream, corresponding to a phys- 
ical length I ~ 6c?3 pc. Assuming most photons in the re- 

5 Although there is recent revision to this formula, e.g. 
Beck & Krausc (2005), we note that the revised version does not 
work for flat spectrum objects in which a > —0.5. 



gion have frequencies v < 10 16 Hz (w 0.04 keV, nearly the 
detection threshold of Chandra), the synchrotron cooling 
time scale 



' syn 



> 2000 



f B \ 



\20nGJ U0 16 Hz 



yr 
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suggests a flow speed 
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Wfiow,R 



< 3000 * {m)' 



10 16 Hz 



kms 



(4) 

Here we ignored adiabatic cooling, since the nebula cross- 
section shows no widening beyond the central bulge. 
Moreover, if the majority of particles radiate below the 
break frequency of 10 13 Hz, the flow speed will then be as 
low as ~ 100 km s -1 . We emphasize that these estimates 
are highly uncertain because the exact spectrum is un- 
known. Nonetheless, a comparison to the X-ray results 
suggests a significant deceleration of the flow when mov- 
ing downstream, with the speed dropping by 1-2 orders 
of magnitude. 

4.3. Magnetic Field Geometry and Physical 
Interpretation 

The most remarkable feature of G319.9-0.7 is the dis- 
tribution of the polarization vectors. Since synchrotron 
emission of a PWN is optically thin and Faraday thin 
(i.e. depth depolarization is negligible) at our observed 
wavelengths, the polarization maps in Figure |4] repre- 
sent the magnetic field of the nebula projected on the 
sky plane. This implies a helical field structure in the 
north, with the axis of symmetry parallel to the nebu- 
lar axis, and hence to the pulsar's inferred direction of 
motion. This is the first time an azimuthal field geome- 
try has been seen in a bow-shock nebula. Axisymmetric 
toroidal fields about pulsar rotation axes have been ob- 
serve d in young PW Ne within SNRs (e.g. iDodson et al.l 
2003; Hester 200|), and have b een success fully repro- 
duced by MHD simulations (see, iVolpi et al.1 12008L and 
references therein). As an analogy to these systems, the 
field configuration of G319.9-0.7 could be understood if 
the pulsar motion aligns wi th its spin axis. Su ch an 
p alignment is not uncommon (|Johnston et al.l [20051 ). and 
'could be the resu lt of a momentum kick at the birth of a 
neutr on star (see. lSpruit fc Phinnevlll998tlNg fc Romanil 
120071 ). If this is the case for PSR J1509-5850, the jet- 
like polar outflow would contribute to the X-ray tail. The 
MHD instabilities (e.g. kink and sausage instabilities) in 
the flow could then help explain t he nonuniform surface 
brigh tness in the tail downstream (|Kargaltsev fc Pavlovl 
l2008f ). Future polarization measurements of the pulsar 
profile may reveal the projected orientation of its spin 
axis, directly confirming the above picture. However, 
since PSR J1509-5850 has a relatively lo w degree of po- 
larization ()Weltevrede k, Johnston! [2008). this may re- 
quire the next generation of radio telescopes. 

Further support for the spin-velocity alignment is 
given by analyt ic stud ies of pulsar wind bow-shocks. 
iRomanova et al.l ((2005) showed that an aligned pulsar 
can wind up its magnetic field to form an azimuthal field 
in the magnetotail behind, possibly resembling what we 
found in G319.9-0.7. Assuming energy equipartition be- 
tween the field and the particles, these authors derived 
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an analytic expression for the magnetotail radius r m (z) 
in terms of the distan ce z from the star (Eq. 34-35 in 
iRomanova et al.l 120051. Using the bow-shock standoff 
distance of 0'.'5 (|K08l ). their model predicts a too nar- 
row nebular morphology. If the model parameters are 
allowed to vary, we are able to obtain a better fit, which 
is overplotted in Figure [H although it requires a very 
large standoff distance z s h « r m (0) = 2.0 x 10 17 d3cm 
corresponding to 4 ", eve n at the lower limit 5000 km s 
of the flow speed (|K08t ). The fit also suggests a large 
deceleration in the flow of length scale L = 3 pc, much 
shorter than the total distance traveled by the pulsar, 
as proposed in their model. This could explain why the 
observed nebular morphology is much wider than the the- 
oretical predictions. 

Beyond the central bulge, the magnetic field of G319.9- 
0.7 switches direction abruptly and runs parallel to the 
ne bular axis, similar to that of the Mouse PWN reported 
bv lYusef-Zadeh fe Gaenslen pfiOST ). The change in field 
geometry seems difficult to understand. One possible 
scenario could be a two-component field structure con- 
sisting of a poloidal core surro unded by a toroidal fi eld, 
such as the one suggested by iKothes et al.l (|2008l ) for 
DA 49ifl If the toroidal field decays faster downstream, 
this would leave the poloidal core component, resulting 
in a switch in the overall field orientation. In this pic- 
ture, the projection of the two field components leads to 
a lower degree of polarization in the PWN interior than 
around the edge as observed. Near the transition region 
where the two field strengths are comparable, we expect 
the total polarization to cancel out. However, such a 
complete depolarizing region is not observed in Figure O 
Indeed, Figure [4] shows that the polarization vectors near 
the bulge have intermediate orientation, suggesting phys- 
ical change of the field direction rather than projection 
effect. Therefore, we conclude that this scenario is less 
likely. 

Hinted by the flow deceleration and the relatively 
lower degree of polarization at the bulge, an alterna- 
tive physical scenario is that the plasma flow might be- 
come unstable downstream, then drives an expanding 
bubble into the ISM, s imilar to the model suggested by 
Ivan Kerkwiik fe Ingld ([2008) to explain the Ha observa- 
tions of the Guitar Nebula, the bow-sh ock PWN powered 
by PSR B2224+65 (jCordes et al.lll993l ). If we decompose 
the magnetic field into random and ordered components, 
then the relative strength q 2 — ^a't, between the ran- 
dom field b and the ordered field perpendicular to the 
line of sight B± can be estimated from the degree of po- 
larization p. For the case of energy equipartition, 



P_ 
Po 



l + 3q 2 



(5) 



where pp = 5,3 ° is the intrinsic degree of polarization 
(|Beck fc Kraiisel 120051) . As indicated in Figure [SJ the 
nebular emission is ~ 40% polarized around the edge, 
implying a highly ordered field of q 2 « 1.2. On the other 

6 We should note that the field structure of DA 495 is inferred 
from the RM distribution of the PWN. However, internal Faraday 
rotation is generally not expected for non-thermal plasmas such as 
pulsar winds. 



hand, the lower degree of polarization p ~ 25% in the in- 
terior gives q 2 ~ 3.2. This indicates a substantial random 
field, suggesting some hints of turbulence, which could be 
due to flow instability. If this is the case, then the ac- 
tual value of q 2 would be even h igher due to the i ntrinsic 
anisotropy of MHD turbulence (jBeck et al.ll2003l ). 

Assumin g continuous energy injectio n and adiabatic 
expansion , Ivan Kerkwiik fc Ingli (12001 give the angular 
radius b of the bulge in terms of the instability angular 
scale A using the Sedov- Taylor solution: 



,2/5 




(0d - A) 



3/5 



(6) 



where / b is a function of order unity, v* is the pul- 
sar velocity, 0d is the angular separation between the 
bubble center and the pulsar, and #0 is the angular 
size of the bow-shock stan doff distance (see Eq. 4 in 
van Kerkwiik fc I ngle 2008). Our radio maps in Figure [TJ 
indicate b = 80" and d = 280", together with O = 0"5 
and w 400 km s" 1 from IK08I : these suggest A « 170", 
implying a relatively slow growth rate of the instabil- 
ity. The ratio A/0o ~ 340 is higher than the value 80 
inferred from the Guitar Nebula, but still within a few 
times of the instability growth leng t h scal e A g < 60 9q as 
discussed bv lvan Kerkwiik fc Ing ic (2008). A more real- 
istic estimate with the magnetic field taken into account 
could lead to a lower expansion rate, and hence a slightly 
smaller A. The estimate above predicts an expansion rate 



2d» b 



Mb 



5 d d - A 



= 8 mas yr 



(7) 



for the bulge, too small to be detected by the ATCA over 
a reasonable time span, and the emission is too faint for 
VLBI measurements. 

For the orientation switch of the polarization vectors 
downstream, we speculate that the expanding bubble 
may channel particles into the low-pressure cavity evac- 
uated by the pulsar wind along the path, thus changing 
the field configuration. Finally, we should point out that 
this picture fails to explain how the polarization struc- 
ture is preserved during the bubble expansion. Further 
simulation work is thus necessary to complete the study. 

4.4. Multiwavelength Comparison 

Figure [3] show that the X-ray tail is brightest near the 
pulsar while t he rad io emission is extremely faint in the 
same region. IK08I extracted the X-ray counts from a 
14" x 24" ellipse at the head of the PWN (excluding the 
pulsar), and obtained an X-ray power-law spectrum with 
index a = —0.8 ± 0.3. If the absorption column density 
is allowed to vary, they found a = —0.27 ±0.25, identical 
to the radio spectral index of the overall PWN. Using 
the same extraction region, we obtained flux densities of 
0.5 ± 0.3 mJy and 0.6 ± 0.2 mJy at 3 and 6 cm, respec- 
tively, from the radio intensity maps in Figure [1] Given 
the X-ray spectrum, the faintness in radio emission is 
not unexpected. If we extrapolate the two different X- 
ray power-laws to radio wavelengths without any spec- 
tral breaks, we expect flux densities of 0.5 /iJy-14mJy at 
3 cm and 0.6 /xJy-20 mJy at 6 cm, consistent with the ob- 
served values, albeit not very constraining. To determine 
if there is any intrinsic break in the injection spectrum, 
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further studies will require deeper radio and X-ray ob- 
servations, as well as flux measurements at wavelengths 
in between (e.g. infrared). 

Moving downstream from the pulsar, while the X-ray 
tail fades due to synchrotro n bu rn-off, the radio emissiv- 
ity gradually increases. As IK08I point out, this could be 
caused by the deceleration of the flow, and similar spa- 
tial radio-X-ray anti-correlations have been observed in 
other PWNe with mildly relativistic jet-like polar out- 
flows, e.g. the ones powered by PSRs B1509-58, B1706- 
44 and Vela dGaensler et aljl200£ iRomani et all 120051 : 
iKargaltsev fc Pavlovi^OOl TFor the case of G319.9-0.7, 
the bulge is ~ 15 times brighter than the tip from the 
6 cm linear profile (Figure [3]) and 10 times wider (Fig- 
ure [1]). Therefore, the volume emissivity at the bulge is 
15% of that at the tip. As discussed above, the flow speed 
drops by 10-100 times between the two regions. After ac- 
counted for the geometry, this corresponds to a factor of 
1-10 decrease in particle number density n downstream. 
Due to the strong dependence of synchrotron emissivity 
on the field strength (L/V oc nB 1 1 2 for energy equipar- 
tition, or oc n B 2 if the particle distribution remains un- 
changed), our results indicate a minimal field decay be- 
tween the tip and the bulge, at most by a factor of 2. 
While this seems incompatible with the field strengths 
estimated in §4.21 it argues that energy equipartition is 
not achieved near the pulsar and the actual magnetic 
field could be close to 20 fiG at the tip of the PWN. 

4.5. Comparison with Other Radio Bow-shock PWNe 

For a handful of bow-shocks with detailed radio stud- 
ies, their radio emissions generally peak at the pulsar 
position and show good spatial correlations with the X- 
ray counterparts, e.g. the Duck, the Mouse and IC 443 
(|Kaspi et al.l[200lt iGaensler et alJl200ll2006ft . The clear 
radio-X-ray anti-correlation of G319.9-0.7 likely reflects 
a different flow condition for the source compared to 
these other systems. Additionally, an azimuthal -B-field 
geometry has never previously been found in bow-shocks. 
In the following discussion we will focus on the Mouse 
PWN, since it has high resolution polarimetric observa- 
tions that allows a direct comparison with our study. 

The magnetic field of the Mouse wraps around the 
bow-shock at the apex, it then switches orientation 
abruptly behind the pul sar and runs parallel to the neb - 
ular axis downstream ([Yusef-Zadeh fc Gaenslerl [2005). 
Since the radio emission is very faint at the tip of G319.9- 
0.7, it is unclear if a similar field direction switch may 
exist near the pulsar. The azimuthal field geometry in 
G319.9-0.7 is not found in the Mouse. This could reflect 
different spin-velocity alignments of the central pulsars 
in the two systems. While we suggested that G319.9-0.7 
could be an aligned case, we may expect a completely dif- 
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ferent field geometry if the misalignment angle is large. 
Future MHD simulations with anisotropic winds and pul- 
sar orientations taken into consideration could help verify 
this picture. 

More intriguing is the field structure at large distances 
from both pulsars. Figure [4] shows that the field orienta- 
tion of G319.9-0.7 exhibits a 90° switch near the central 
bulge, at ~ 4.5 e?3 pc from the pulsar. Similar behavior is 
also found in the Mouse: the polarization vectors change 
direction abruptly at several places along the nebular 
axis, with the most prominent one at ~ 5 d§ pc down- 
stream from th e pulsar, and the behavior appears to be 
quasi-periodic (Yusef-Zadeh & Bally [19871) . This could 
be due to MHD instabilities in the flows, but the detailed 
mechanism is not well understood. More examples are 
necessary to determine whether this is a common feature 
among bow-shocks with long trails. 

5. CONCLUSIONS 

We have presented a detailed radio polarization study 
of G319.9-0.7 at 3 and 6 cm. Our results confirm 
that this source is a bow-shock PWN associated with 
PSR J1509-5850, and reveal its central bulge morphol- 
ogy, which could represent an adiabatically expanding 
bubble drove by the flow instability. We found an 
equipartition field strength of ~ 20 /iG and a flow speed 
< 3000 km s -1 in the radio nebula, indicating a substan- 
tial deceleration of the flow as compared to the X-ray- 
emitting region near the pulsar. The polarization mea- 
surements reveal a remarkable azimuthal field geometry 
in the northern half of the PWN, which qualitatively 
agrees with the theoretical prediction for an axisymmet- 
ric bow-shock. As compared to the well-studied case of 
the Mouse PWN, G3 19. 9-0. 7 shows very different field 
structure and surface brightness distribution, illustrating 
the diversity in the physical properties of bow shocks. 

To conclude, our study demonstrates that high resolu- 
tion radio polarimetry can provide a powerful diagnostic 
tool for probing the physical conditions of a PWN. Future 
works on expanding the sample of polarimetric observa- 
tions on bow-shock nebulae are essential to understand 
the general properties of the population. 
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Fig. 1. — Radio continuum images of G319.9-0.7 at (a) 6 cm and (b) 3 cm. The position of PSR J1509— 5850 is marked by the cross. The 
grey scales are linear, ranging from -0.7 mjy beam -1 to +6mjybcam~ 1 in the 6cm map, and from -0.8 mjy beam -1 to +5mjybeam -1 
in the 3 cm map. The beam size of FWHM 17" X 15" is shown at the lower left of each panel. 
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Fig. 2. — Exposure-corrected Chandra X-ray image of the field around G319.9— 0.7 in the 0.5— 7keV energy band, smoothed to 6" 
resolution and overlaid with the 6 cm radio intensity contours from Figure (TJ a) at 0.5, 2, and 4.5 mjy beam -1 . T he cr oss marks the 
position of PSR J1509— 5850 and the dotted line shows the best-fit theoretical model given by Romanova ct al. (2005) (see 814.3^ . 
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Fig. 3. — Linear brightness profiles of G3f9.9-0.7 in radio and X-rays, extracted from 3' wide rectangular regions along the nebular 
axis. The dotted and solid lines represent the radio profiles obtained from the 3 and 6 cm intensity maps in Figure [T] and the histogram is 
extracted from the exposure-corrected Chandra X-rays image binned to 25". 
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Fig. 4. — Intrinsic orientation of polarization B- vectors of G319.9-0.7 at (a) 6 cm and (b) 3 cm, after correction for the foreground Faraday 
rotation, overlaid with the total intensity contours at the corresponding wavelengths. The 6 cm map is clipped at a signal-to-noise ratio 
(S/N) of 5 for the polarized intensity, and the 3cm map is clipped at a S/N of 3. The contours are at levels of 0.5, 2, and 4.5 mjy beam -1 
for the 6 cm map and 0.6, 2, and 4mjybeam~ 1 for the 3 cm map. The vector lengths are proportional to the polarized intensity, with the 
scale bar at the lower left representing 1 mjy beam - 1 . The position of PSR J1509-5850 is marked by the cross, and the beam size is shown 
at the lower left of each panel. 
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Fig. 5. — Fractional linear polarization map of G319.9-0.7 at 6cm, clipped at a signal-to-noise ratio of 5 for the polarized intensity, 
overlaid with the total intensity contours at levels of 0.5, 2, and 4.5 mjy beam -1 . The grey scales are linear, ranging from 10 to 52%. The 
position of PSR J1509-5850 is marked by the cross and the beam size is shown at the lower left. 
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Fig. 6. — Map of rotation measure (RM) towards G319.9-0.7 obtained from the 20cm SGPS data, overlaid with the 6cm intensity 
contours at 0.5, 2, and 4.5 mjy beam -1 . The filled and open squares indicate positive and negative RM values, respectively, with the 
linear dimensions of each box representing the RM magnitude, ranging from -100 to +70radm -2 . The cross marks the position of 
PSR J1509-5850, and the beam size of the 20 cm observation (FWHM of 72" X 63") is shown at the lower left. 
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Fig. 7. — (a) Radio continuum spectrum of G319.9— 0.7, after subtracting the compact source in the south. Due to limited spatial 
sampling, the flux at 3 cm is likely underestimated (see text). Fitting a power-law spectrum to the other three data points suggests a 
spectral index a = —0.26 ± 0.04. (b) Broadband spectrum of the nebula by comparing to the X-ray results from [K08l . The dotted lines 
indicate the uncertainty in the X-ray power-law spectrum. 
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